Theory of Half-metallic Ferrimagnetism in Double Perovskites 
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We present a comprehensive tiieory of tiie temperature- and disorder-dependence of lialf-metallic 
ferrimagnetism in tlie double perovskite Sr2FeMo06 (SFMO) with Tc above room temperature. We 
show that the magnetization M{T) and conduction electron polarization P{T) are both propor- 
tional to the magnetization MsiT) of localized Fe spins. We derive and validate an effective spin 
Hamiltonian, amenable to large-scale three-dimensional simulations. We show how M{T) and Tc 
are affected by disorder, ubiquitous in these materials. We suggest a way to enhance T^ in SFMO 
without sacrificing polarization. 



Double perovskites (DPs) A2BB'06 are an important 
family of complex oxides, derived from the simple ABO3 
perovskite structure by a three-dimensional (3D) checker- 
board ordering of B and B' ions. One of the best stud- 
ied examples is Sr2FeMo06 (SFMO), a half-metaUic fer- 
rimagnet with Tc ~ 420K, well above room tempera- 
ture prtil- Clearly SFMO, and other DPs, can have 
enormous technological impact if their stoichiometry and 
ordering can be controlled. 

From a theoretical point of view, we argue that DPs 
are simple systems for understanding metallic ferromag- 
netism, despite their apparent complexity. First, in con- 
trast to iron, there is a clear separation of the local- 
ized (B) and itinerant degrees of freedom (coming from 
B') in the DPs. Second, in contrast to the mangan- 
ites, DPs have neither Jahn- Teller distortions nor com- 
peting superexchange, given the large distance between 
B-sites. Third, in contrast to dilute magnetic semicon- 
ductors, disorder is not an essential aspect of the theoret- 
ical problem. Important early theoretical work on half- 
metallic DPs includes r=0 electronic structure calcula- 
tions 3] , and model Hamiltonians analyzed using various 
mean-field theories [J-Isl and two-dimensional (2D) sim- 
ulations |7|. 

In this Letter, we present a comprehensive theory 
that gives insight into the temperature- and disorder- 
dependence of the magnetic properties of metallic DPs. 
We make detailed comparisons with and predictions for 
SFMO. Our main results are: 

(1) We show that both the total magnetization M{T) 
and the conduction electron polarization P{T) at Ef are 
proportional to the magnetization Ms{T) of localized Fe 
spins. This result is significant because, while M{T) is 
easy to measure, it is P{T) that is of crucial importance 
for spintronic applications. 

(2) Our main theoretical advance is the derivation and 
validation of an effective classical spin Hamiltonian HoS 
[see eq. ^] for DPs, which differs from both the Heisen- 
berg and Anderson-Hasegawa models [8| . We show that 
Hcs describes the full T-dependence of the magnetiza- 
tion Ms{T), and hence that of M{T) and P{T). 

(3) We present the results of simulations of iJcff on large 



3D lattices, including disorder effects, thus going beyond 
all previous theoretical calculations on SFMO. 

(4) We compute M{T) and Tc, using microscopic band- 
structure parameters as input, and see how these are af- 
fected by deviations from stoichiometry and by anti-site 
(AS) disorder, ubiquitous in real materials. Ours is the 
first theory to show that Tc is insensitive to AS disorder, 
in excellent agreement with experiments, even though 
Af (0) is suppressed. 

(5) We conclude with a novel proposal to enhance Tc of 
SFMO without sacrificing polarization, using a combina- 
tion of disorder and doping. 

Model Hamiltonian: For large Hund's coupling Jh, 
the Fe^+ (3d^) site has a S=5/2 "core spin" or local 
moment. The Mo^+ (4^^) contributes a t2g electron 
which hybridizes via O with the Fe t2g states. Sym- 
metry implies that da/3 electrons delocalize only in the 
{a, /3)-plane l(j|. Thus the motion of electrons in the 3D 
system decouples into three 2D planes. The "double ex- 
change" Hamiltonian [5|, LZj describing itinerant electrons 
interacting with core spins is: 

if = -f ^ {eicrdl^Cja + h.c.) 

-t' E 4c,'. + ^E<^»^ (1) 

Here dia (cia) are fermion operators on the Fe (Mo) 
sites with spin a. At the Fe sites i, we choose lo- 
cal axes of quantization along Si and Pauli exclusion 
prohibits an f electron. For all the Mo sites j, we 
choose the same (global) axis of quantization. The ori- 
entation {9i,ipi) of the classical spins Si then affects 
the Mo-Fe hopping via Ci^ = — sin(6'i/2) exp{i(pi/2) and 
£-4 = cos(6'j/2)exp(i(^i/2). 

The parameters in H are the hopping amplitudes t, 
between nearest neighbor (Fe-Mo) sites, and t', between 
two Mo sites, and the charge transfer energy A between 
Fe t2g4, and Mo t2g states; see Fig. [Ija). The results 
are independent [9| of the sign of t. For now, we choose 
t = 1 as our unit of energy. Symmetry dictates i' > 0. 
We choose t' /t = 0.1 and A/t — 2.5, using realistic band 
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FIG. 1: (a) Schematic showing energy levels at transition 
metal sites in two unit cells (formula units) of SFMO. The Fe 
sites have localized S — 5/2 core spins, treated as classical 
vectors with orientation {9,ip). The parameters t,t' and A 
of the Hamiltonian ([1} , governing the dynamics of the itiner- 
ant electrons in t2g orbitals, are also shown, (b) Calculated 
electronic structure -E(k) and the spin-resolved DOS Nf{E) 
and Ni{E) in the ferromagnetic ground state with all core 
spins up. Note the half-metallic ground state in SFMO with 
conduction electrons polarized opposite to core spins. 



parameters for SFMO as input [3|. We will show below 
that t sets the scale for the magnetic Tc, and t = 0.27 
eV, consistent with ref. [3|, leads to the experimental 
Tc = 420 K of SFMO. 

We treat the quantum mechanics of "fast" itinerant 
electrons using exact diagonalization (ED) in the back- 
ground of "slow" S=5/2 spins, for which we use a T 7^ 
classical Monte Carlo (MC) simulation. This separation 
of time scales is justified below. Using standard ED- 
MC techniques, supplemented by T—0 variational calcu- 
lations, we obtain the results shown in Figs. [T] and [5] 
These are limited to 2D systems for computational rea- 
sons. 

We show in Fig. [IJb) the band structure and the spin- 
resolved density of states (DOS) N„{E) for the ferri- 
magnetic ground state for the band-filling correspond- 
ing to SFMO with n = 0.33 electrons per unit cell per 
plane. Here the conduction electrons, with magnetiza- 
tion Mei{T = 0) = l/is (per unit cell), are polarized 
opposite to the Fe spins, with core spin magnetization 
-^S'(O) = 5/iB. The half-metallic ground state, with 
iV^(O) = and A'^(O) 7^ 0, has a net magnetization 
M{0) = Ms{0) ~ M,i{0) ^ 4nB- 

We see from Fig.UJa) that Ms{T), Mei{T), and hence 
the total magnetization M{T) = Ms{T) - M^iiT), all 
have essentially the same T-dependence. Another quan- 



tity of great interest is the conduction electron polariza- 
tion P{T) at Ef which determines the tunneling mag- 
netoresistance. This is defined by P{T) = [N^(0) — 
iVt(0)]/[A^i(0) + iVt(O)]. with ^ = in the DOS is mea- 
sured from Ef. From Fig. [IJb), we find that P{T) also 
follows the T-dependent Ms{T). 

Effective Spin Hamiltonian: The results of 
Fig. [2ja) and (b) imply that if we had a reliable theory 
for the core spin magnetization Ms{T), we could under- 
stand all the magnetic properties of SFMO, including the 
polarization P(T). With this motivation, we derive an 
effective Hamiltonian Hcs for the core spins by gener- 
alizing the two-site Anderson-Hasegawa [8] analysis for 
manganites to double perovskites. 

To derive i?off, we find the exact solution of ([T|) for 
two unit cells. The Hilbert space, for a given t2g orbital, 
has three states per unit cell: Fe t2gi and Mo t2g,^,i. 
The smallest accessible filling i s n = 0.5 (per plane), i.e.. 



one electron in two unit cells [ll|. We analytically find 
the lowest eigenvalue as a function of the angle {6i — 
6j) between spins. Working in two different geometries, 
we find [12|| the nearest-neighbor (Ji) and next-nearest- 
neighbor (J2) interaction energies; see inset in Fig. [5Jc). 
Expressing these in terms of Si • Sj, where each Si is a 
unit vector, we obtain the effective Hamiltonian 



^cff — —Ji 



(iJ) 



Fi (S, • S,) - J2 5] F2 (S, • S,) 



(2) 



where the functions Fi{x) ~ By 2 -f- y/2 + 2x and 
F2{x) = {5 + V5)^/6 + 2V3 + 2x. Our two- unit cell anal- 
ysis gives explicit expressions [12| for Ji and J2, both of 
which are ferromagnetic with their scale set by the ki- 
netic energy t of delocalization. We emphasize that the 
double square-root form of HcS is quite different from 
the (single square-root) Anderson-Hasegawa model. 

Next we need to understand how we can use iJoff going 
beyond the simple two-unit cell derivation. Specifically: 
(i) How can we relate Ji, J2 to t,t',A and the filling 
n? (ii) To what extent does HcS capture the essential 
physics of the full Hamiltonian H7 

The dependence of Ji and J2 on microscopic parame- 
ters can be obtained by matching the spin-wave (SW) 
spectra of Hcs and H. This comparison is shown in 
Fig.HJc) along certain symmetry directions. We find SW 
dispersion for H using ED to compute the energy of elec- 
trons moving in a "frozen" spin- wave background. The 
SW analysis for Hcs is straightforward, since for small 
deviations from the FM ground state, Hcs reduces to a 
nearest and next nearest neighbor FM Heisenberg model 
i^Hcis- The low energy scale of O.li for spin dynamics (see 
Fig-Hlc)) justifies a posteriori our assumption of "slow" 
spins and "fast" electrons, whose bandwidth is of order t 
(see Fig. HUb) ) . The same separation of energy scales also 
justifies the use of T- independent exchange couplings Ji 
and J2 for all T < Tc, as we discuss next. 
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FIG. 2: (a) Magnetization of the core spins Ms{T) and of 
the conduction electrons Mei{T) calculated using ED-MC 
method, (b) Conduction electron polarization (at Ef) P{T) 
and normalized Ms{T). We see that all magnetic properties 
to be proportional to Ms(T). (c) Spin wave dispersion of full 
Hamiltonian H , obtained by exact diagonalization, compared 
with that of the effective Hamiltonian HeS- Inset: Fe-Mo lat- 
tice showing nearest-neighbor (Ji) and next-nearest-neighbor 
(J2) interactions of -ffcft- (d) The normalized magnetization 
M{T) for three Hamiltonians: H, Hc« and the Heisenberg 
model. All results are obtained on 8 x 8 systems with error 
bars no larger than the symbol size. 



To validate the effective classical model H^s we show in 
Fig. [2jd) that it reproduces the magnetization M{T) of 
the full Hamiltonian H over the entire range of tempera- 
tures. We note that the corresponding Heisenberg model 
^Hcis gives quite different results, except in the T — >■ 
limit with small spin deviations. In other words, the T- 
dependent magnetization of DP's cannot be modeled by 
a Heisenberg Hamiltonian. But the (highly non-linear) 
iJeff that we have derived gives an excellent description 
of the ED-MC result for the full H. 

The classical Hes can be easily simulated on large 3D 
lattices, unlike the full iJ, and the results are shown in 
Fig. [3 We note the linear drop in M{T) at low T, due to 
classical spin- waves, followed by a rapid suppression of M 
at the phase transition I13| . We estimate Tc in the infinite 
volume limit using finite size scaling; see Fig. [3]Jb). For 
t' It = 0.1 and A/i = 2.5, we find T^ = 0.14i. Comparing 
this to Tc = 420K for pure SFMO, we obtain t = 0.27 
eV, consistent with ref. [3J. 

Disorder: iJeff permits us to model various kinds of 
disorder and deviations from stoichiometry [12|. (i) Ex- 
cess Fe is modeled with an extra spin (at a Mo site) 
that interacts with its neighboring spins with a large an- 
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FIG. 3: (a) Magnetization M{T) from 3D simulations of HaB 
on 16"^ systems. The infinite system Tc is obtained from the 
scaling analysis of panel (b). (b) Finite size scaling of data for 
L'-^ systems. Here e = \T - Tc\/Tc and /3 = 0.36 and ly = 0.70 
are the well-known 3D 0(3) critical exponents. 
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FIG. 4: (a) M{T), normalized by the M(0) for the disorder- 
free system, for different degrees of anti-site (AS) disorder 5 
(see text), (b) Comparison of theoretical results for Tc and 
r=0 magnetization A/(0) with experiments [13|. M{0) drops 
like 2S, while Tc is insensitive to AS disorder. 



tiferromagnetic (AF) superexchange S{S + 1) Jaf — 34 
mev [ij]. (ii) Excess Mo is modeled by removing an Fe 
spin from the lattice. Both (i) and (ii) also require reeval- 
uation of Ji and J2 due to change in carrier density n. 
(iii) Here we focus on anti-site (AS) disorder, the most 
common form of disorder in DPs, with Fe and Mo inter- 
changed and no change in n. 

We see from Fig. U that AS disorder systematically re- 
duces M(0) without affecting T^ in excellent agreement 
with experiments [la, II4I • We quantify AS disorder using 
S the fraction of Fe atoms that are on the Mo sublattice. 
The observed magnetization M(0)[1 — 25] arises from the 
loss of two moments for each AS defect. One from the 
moment lost at the Mo (on the Fe sublattice) and the 
other from the Fe (on the Mo sublattice) antifcrromag- 
netically coupled to its neighbors via Jaf- 

There are two opposite effects of AS disorder on Tc that 
appear to balance each other. The strong Fe-Fe superex- 
change Jaf pins the spins surrounding the Fe-defect, and 
makes the magnetic order more robust against thermal 
fluctuations. On the other hand, the Mo-defect leads to 
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FIG. 5: (a) Core spin magnetization Ms{T) for 
Laa;Sr2-a;Fei+yMoi_j,06 for uncompensated {x,y) = (0,0.25) 
and compensated {x,y) = (0.75,0.25) Fe-rich systems, (b ) 
The uncompensated Tc{y) is compared with experiments [2y|. 
Our prediction for compensated {x = 3j/) system shows large 
enhancement in Tc- 



broken Ji , J2 bonds which weaken the magnetism. The 
net effect is a T^ insensitive to 5 for moderate levels of 
AS disorder, which is exactly what experiments observe. 
We note that while the loss of M(0) with AS disorder has 
been explained earlier [16l ll8[ ll9| , ours is the first theory 
to correctly account for Tc] previous theories either found 
a drop [16[ in Tc or an increase [5(. 

Raising Tc- We conclude with a proposal to raise 
the ferromagnetic Tc without sacrificing the conduction 
electron polarization P; see Fig. [5] In brief, it involves 
adding excess Fe and compensating for the loss of mobile 
carriers by La substitution on the Sr site. 

To see how this works, let us first consider excess Fe 
without any compensation: Sr2Fei+j^Moi_j,06. In this 
case, Ms{Q) decreases like y due to the AF alignment of 
excess spins. Tc also decreases with y because the loss 
of carriers n = (1 — 32/)/3 dominates over the enhanced 
pinning of moments by Jaf in the vicinity of defects . For 
the uncompensated case, the calculated Tc{y) in Fig.ISfb) 
is in reasonable agreement with experiments [2(1] . 

We can compensate for the carriers by A-site substi- 
tution: La2;Sr2._xFei+yMoi_j^06 with electron density 
n = {\ + X — 3y)/3 per plane. Choosing a; = 3j/ counters 
the doping-dependent drop in Tc that dominated above 
and we find that Tc can be significantly enhanced over 
that of pure SFMO due to just the local pinning of mo- 
ments at excess Fe sites. This increase in Tc goes hand- 
in-hand with an unchanged Ms(0); see Fig. [5l^a). 

We also note that the polarization P(0) remains 100%. 
When the Fe spins cant at finite T, electrons depolar- 
ize by mixing of up and down states, however, no such 
processes are permitted at T=0, even in a disordered Fe- 
rich system. We have checked using the full Hamiltonian 
dl]) that, even though the strict proportionality between 
Ms{T) and P{T) is not observed in the presence of dis- 
order, the high temperature P is still enhanced over the 
clean system due to the large increase in Tc- 

An alternative way to enhance Tc is to add mobile elec- 
trons using La doping, which, however, has been shown 



to lead to considerable increase in AS disorder [15|. In 
principle, compensated doping as proposed here should 
introduce less AS disorder [21 1. 



In conclusion, while we have focused here on SFMO, 
our theory provides a general framework for understand- 
ing half metallic ferrimagnetism in DPs. Interesting 
directions for future work include A-site substitution. 
Coulomb correlations on B', which may become increas- 
ingly important for larger carrier concentrations, and 
spin-orbit coupling on B' for 5d elements. 
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